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Abstract. We present model light curves, Power Spectral Densities (PSD) and time 
lags of accreting Black Hole Candidates (BHC) based on a recent model of these 
sources. According to our model the observed variability is due, to a large extent, 
to the stochastic nature of Comptonization, the basic process responsible also for the 

^ ' emission of high energy radiation. Our additional assumption is that the Comptoniza- 

tion process takes place in an extended but non-uniform atmosphere around the com- 
pact object. Our model reproduces the observed light curves well, in that it provides a 
good fit to the PSDs and the overall light curve morphology, indicating, in accordance 
with observation, that most of the power resides at time scales <; a few seconds while 
at the same time one can distinguish shots of a few msec in duration. We suggest that 
0^ ' refinement of this type of model along with spectral and phase lag information can be 

(— I used to probe the structure of this class of high energy sources. 
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The study of the physics of accretion onto compact objects (neutron stars and 
^ black holes) involves length scales much too small to be resolved by the current or 

the forseable future technology. As such, this study is conducted mainly by the 
theoretical interpretation of spectral and temporal observations of these systems. 
In this respect, spectral analysis presents the first line of attack in unfolding their 
physical properties. For the specific systems and in particular Black Hole Candidate 
(BHC) sources, a multitude of observations have indicated that their energy spectra 
can be fitted very well by Comptonization of soft photons by hot electrons. It is 
thus generally agreed that Comptonization is the process by which the high energy 
(> 2 — 100 keV) spectra of these sources are formed. Thus, while the issue of the 
detailed dynamics of accretion onto the compact object is still not resolved, the 
assumption of the presence of a thermal distribution of hot electrons has proven 
sufficient to produce models which successfully fit the spectra of the emerging high 
energy radiation. Spectral fitting, has been hence used as a probe of the dynamics 
of accretion onto the compact object. 
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INTRODUCTION 



However, the Comptonization spectra, which depend roughly on the product of 
the electron temperature and the photon escape probability from the scattering 
cloud, do not provide in and of themselves a measure of the scale of the system. As 
such they cannot provide any clues about the dynamics of accretion of the hot gas 
onto the compact object, let alone its parameters (density, velocity) as a function 
of radius, necessary to determine these dynamics. One needs, in addition, time 
variability information. 

While simple minded dynamics anticipate most of the variability power to be 
associated with the dynamical time scales of the last few gravitational radii (msec 
for galactic sources), the observed PSD indicate that the variability power resides 
at time scales ^ 1 sec, far removed from the "natural" ones. At the same time, 
however, the observed light curves do exhibit the presence of "flares" with the 
expected (a few msec) rise times (see e.g. [1]), even though the latter make an 
insignificant contribution to the PSD. It is therefore thought that the observed 
PSDs are due to a modulation of the accretion rate onto the compact object, the 
result of a kind of self-organized-criticality for the dynamics of accretion [2] . 

However, the Comptonization process, provides a much more refined probe of 
variability than the PSD: Because, on the average, the energy of the escaping 
photons increases with their residence time in the scattering medium, the hard 
photon light curves lag with respect to those of softer photons by amounts which 
depend on the photon scattering time. Therefore, if the high energy radiation is 
emitted from a compact region near the black hole, the resulting lags At should be 
independent of the Fourier components of the time variation and At ~ msec. 

These lags were sought for and Fourier analyzed in the light curves of Cyg X-l 
[3-5], the X-ray transient source J0422+32 [6] and the galactic source GRS 1758- 
258 [7]. Surprizingly, rather than being constant (~ msec), as expected, their 
magnitude was found to increase linearly with the Fourier period P, from <$.001 
sec at P ~ 0.05 sec to roughly 0.1 — 1 sec for P ~ 10 sec. Their magnitude and 
in particular their P-dependence, exclude the possibility that the observed PSD 
are due to modulation of the accretion rate and strongly indicate scattering in a 
very extended non-uniform medium, whose density profile reflects on the functional 
dependence of the lags on P as suggested by [8-10]. 

As suggested originally in [8] and discussed in more detail in [9-11], the issue 
of time lags can be resolved if the Comptonization takes place in a very extended 
(< 10 11 cm) non-uniform hot plasma with density profile n(r) oc r~~ p , p ~ 1, the 
range in P reflecting the range in radius over which this configuration extends. It 
was further discussed there that the dependence of lags on P relates to the value 
of the index p and it can be used to probe the density structure of the high energy 
emitting region. Herein we discuss and present models on the X-ray light curves 
associated with this specific model of high energy emission from accreting compact 
objects. 
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FIGURE 1. The shot profiles as computed by the Monte-Carlo (dashed line) and their analytical 
fits (solid line), along with the corresponding PSD (dotted line). 



MODELING THE LIGHT CURVES 

The fundamental tenet of our model is that the observed variability is due largely 
to the stochastic nature of Compton scattering. As discussed in [8], the response 
of the extended atmosphere to a 5-function injection of soft photons in its center, 
has a power law (rather than exponential as is the case for a uniform medium) 
time dependence. Our model light curves consist of an incoherent sum of such 
responses, the result of random injection of soft photons at given rate determined 
by the dynamics of the smallest radii, i.e. 

N 

F(t)=J2Qi9(t-Ti) (1) 

1=1 

where g(t) is the response function of the atmosphere, Tj are the arrival times of 
the individual shots and Qi their amplitudes. To avoid introducing information in 
the distribution of Q^s and Tj's we assume that the QiS are constant while the Tj's 
are randomly distributed about a mean value. 

To perform the sum of Eq.(l) we produce an analytic fit to the function g(t) 
whose form has been computed using the Monte-Carlo code developed by one of 
us (XMH). The form used in fitting them is the following: 



9{t) 




if x = t/to < 1; 
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The model and computed response functions are given in Fig. 1. The values of 
the parameters corresponding to the specific figure are to = 0.02, (5 — 10 sec, and 
a ~ 0.50, b—1, 7 = 1.5. The agreement between them is excellent. 

The PSD of the light curve of Fig. 2, computed using Eq. (1), due to the 
independence of the shots, is simply the PSD of an individual shot, also shown in 
Fig. 1. This PSD, besides the segment ex uj~ 2a , reflecting the t a ~ l part of g(t), 
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FIGURE 2. A light curve computed using the prescription given in the text with to = 0.001, 
(3 = 1 sec and a = 0.5. 

exhibits also a flattening at uj < 1/(3 and a cut-off at uj > 10 Hz, due respectively 
to the cut-off and rising parts of the shot. 

In figure 2 we present a light curve given by the above prescription. The 
corrsponding values of the shot parameters are a = 0.5, and t = 0.001, (5 — 1 sec. 
The initial rising segment lasts ~ f$ sec, i.e. the cut-off time of the g(t) function 
and it is only a turn-on transient. Eye inspection reveals the presence of shots 
with a broad range of rising times and durations, down to to, and as long as a few 
seconds, even though only a single type of shot has been used in the model. These 
shots of varying duration are simply the result of the random sum of shots with 
power-law tails. One should note that the specific form of g(t) cannot be discerned 
in their random superposition and one might naively consider the light curve as 
a sum of exponential shots with a distribution of timescales. However, the power 
law structure of the underlying shots is apparent in the PSD of the light curves. 
At the same time, the fact that the time lags, as computed within this simple 
model [10,11], fit the observations in such straigtforward fashion, argues strongly 
in support of the above interpretation. 
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